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Abstract 


Using the framework of the Heavy Quark Effective Theory we have re-analyzed the 
Isgur-Wise function describing semileptonic A, — A, decays in the QCD sum rule approach. 
The slope parameter of the Isgur- Wise function is found to be p? = 1.35 + 0.12, which is 
consistent with an experimental measurement and a lattice calculation. To O(1/m», 1/m.) 
of the heavy quark expansion the integrated A, decay width is used to extract the CKM 
matrix element Va for which we obtain a value of |V| = 0.041+0.004 in excellent agreement 


with the value of |V.4| determined from semileptonic B — D* decays. 
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The study of b — c semileptonic weak decays has been the subject of considerable 
interest in recent years, as a source of information on Væ and as a laboratory for under- 
standing the strong interaction effects and developing nonperturbative QCD methods. A 
considerable amount of work has been carried out in the meson sector, in which the Heavy 
Quark Effective Theory (HQET) [I] was first developed. In the baryon sector a particular 
example is the semileptonic decay A, — A.C. In view of the fact that a recent exper- 
iment from DELPHI [2] shows a discrepancy with the results of previous QCD sum rule 
calculations [3] [4] an updated sum rule analysis for the baryonic A, — A. form factor is 
called for. 

In the heavy quark limit, the hadronic matrix element of the A, — A, transition can be 


simply expressed in terms of a single Isgur-Wise (IW) function defined as follows [5l [6] [7], 
(Ac(v^)|ePb|As(v)) = €(w)ta,(v')Pua,(v) , (1) 


where w = v-v’ is the velocity transfer variable and F is an arbitary gamma matrix. When 
the velocity of the heavy quark changes from v to v' due to the weak decay of the heavy 
quark, the light degrees of freedom undergo a corresponding transition due to their strong 
interactions with the heavy quark. The Isgur-Wise (IW) function (w) is a measure of the 
transition amplitude of the light degrees of freedom. The IW function £(w) is normalized 
to 1 at the zero recoil point w — 1. This value is reduced by a few percent after taking 
into account radiative QCD correction effects [S] D] ! To obtain a theoretical description 
of the whole IW function one must use non-perturbative methods which, at the current 
stage, are beset with large uncertainties. In the decay A, — A./r;, the physical region of 
w lies in the range 1 to 1.43. Usually the IW function is expanded up to the first order in 
w= 

Elw) =1- P(w- 1) + O((w — 1)?) (2) 


The slope parameter f? of the IW function at zero recoil has to be calculated using non- 
perturbative methods. 

The QCD sum rule approach [11] [12] has proven to be a reliable tool to deal with many 
problems in the realm of nonperturbative QCD. It has been used successfully to calculate 
the properties of various hadrons. For instance, besides the light mesons and baryons, 
heavy meson properties were systematically analyzed in the sum rule approach within the 
framework of HQET [8]. In the heavy baryon sector the masses of heavy baryons and 
the IW functions describing their weak transitions were calculated in Refs. [3] [4] [13] [14] 
and [15] [16] [17] [I8] [19], respectively. In Ref. [20], the calculation for the heavy baryons 


began with the full theory and the results of the calculation were expanded in terms of 


! At zero recoil the A, — A. weak transition matrix elements have no 1/mg corrections, and the 1/ m 
corrections are small. 


inverse powers of the heavy quark masses. In the HQET sum rule approach the baryonic 
A, — A. IW function was calculated in BIA], and the slope parameter o? was fitted to lie 
in the range 0.5 — 0.8. However, such low slope values would predict exclusive A, — A, 
semileptonic decay rates dangerously close to the inclusive semileptonic rate [21] [22]. A 
first measurement of the IW function of semileptonic A, — A, transitions has recently 
been reported by the DELPHI Collaboration 2]. The errors on this measurement are 
quite large. Using an exponential parametrisation, they quote a value of p? = 1.59 + 
1.10(stat). When the observed event rates were included in the fit they obtained p? = 
2.03 + 0.46(stat)*U o (syst) [2]. Within the large error bars the experimental slope value 
is compatible with the HQET sum rule results of [3] 4] although the experimental central 
values of [2] are considerably higher than the theoretical sum rule results. 

Theoretically, the above-mentioned HQET sum rule results for the slope parameter 
p? — 0.5 — 0.8 appear to be rather small. Because the number of light quark transitions is 
larger in the heavy baryon case than in the heavy meson case one expects that the slope 
of the baryonic IW function is larger than that of the mesonic IW function. In fact, in 
the large Ne limit, p? will be infinitely large [23]. In the spectator quark model approach 
(6) [71 [24] one finds p? = 20256595 —1/2 which turns into an upper bound p° € 2pfneson — 1/2 
when the interaction between the light quarks is turned on [25] ?. 

Concerning the slope parameter of the mesonic IW function, one finds theoretical values 
of about 1 from sum rule calculations [B]. Experimental numbers for the mesonic slope 
parameter also scatter around 1 261 [27] [28]. Using the spectator quark model estimate one 
thus expects baryonic values of the slope parameter f? in the vicinity of 1.5 or slightly below 
that number. The Skyrme model predicts o? ~ 1.3 [29]. In the infinite momentum frame 
model one has p? = 1.44 and in the relativistic three quark model one finds o? = 1.35 
[I9]. For the baryonic sum rule results, radiative corrections to p? and 1/mg corrections 
to the form factors are not expected to be large enough to solve the discrepancy between 
the large experimental central value for p in [2] and the small QCD sum rule results [3l 4]. 
We therefore concentrate only on the leading order results in our analysis. 

The purpose of this work is to present a new QCD sum rule analysis for the leading- 
order Isgur- Wise function describing the A, — A, transition. In particular, we concentrate 
on the sum rule prediction directly for the slope parameter p?. To start with, we first 
review the sum rule analysis of the two-point Green's function involving two heavy baryon 
currents relevant for the determination of the heavy baryon decay constant and its mass. 
A possible choice of the heavy baryon current with the correct quantum numbers of the 


?The large Ne result in the spectator quark model approach can be worked out to be p? = (Ne — 
1)(peneson — 1/4). Since p?neson > 1/2 according to the sum rule of Bjorken, one then recovers the infinite 
slope result of [23] in the large Ne limit. 


heavy baryon Ag is given by 


J =O Crysqhs , (3) 
where C is the charge conjugation matrix, 7 is an antisymmetric flavor matrix, and h, 
and q are the heavy and light quark fields. Note that there is an alternative choice for the 
heavy baryon current j” which is obtained by the replacement y5 — y5% in Eq.(3). The 
two baryon currents give the same diagonal sum rule for the IW function. 


¿From the correlator 


Pe =a f ate OT G"(@)(0)}10) (4 
one can obtain the two-point sum rule [3l [4] [13] [I4]: 


E l fs E (qq)? -i | (o,GG) 
Oe c M D E E (5) 


where the baryonic decay constant is defined as (0|j’|Ag) = fu. A = mag — mq is 
the binding energy of the Ag-baryon in HQET. T is the Borel parameter. According to 
the duality assumption the higher resonance and continuum contributions to the Green's 
function are approximated by the perturbative contribution above a given threshold s. 
Note that the factor 2 on the l.h.s. of Eq. (5) is missing in the calculation of Ref. [I5]. The 
missing factor 2 is, however, of no relevance for the calculation of À and the IW function. 


To obtain the IW function, one considers the three-point correlator 
Ele eu) =? f drido 0175" (x) hur(OThe(0)*(@2)10), —— (6) 


where e = v- k and € = v' - k'. In the hadronic language the three-point function Eq. (6) 
can be expressed in terms of insertions of hadronic states. The lowest contribution involves 
the IW function, 

PEw) lef .l- 


&(e, €, w) = —— M———. po + resonances . (7) 


(A—e(A—e) 2 


On the other hand, the three-point correlator € can be calculated by using the Operator 
Product Expansion. The perturbative contribution can be written in terms of the double 


dispersion relation, 


pert / - = dé a dë e(€, E, w) 8 
(ec) = f af eem (8) 


with the following spectral density, 


PE: 1 
(29644 (w2 = 1)5/2 


(e?+€7—Qwee')*tr(r7") O(c)O(e’)O(Qwee’—e?—-€”). 


(9) 


1+ y. 1 
e(e, €, w) ve eee 
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The condensate contributions will be included as a series of vacuum expectation values of 
operators ordered by their dimension. 

A systematic uncertainty of QCD sum rule calculations lies in the treatment of higher 
state contributions to the hadronic side of =. Generally the quark-hadron duality assump- 
tion is adopted, which simulates the higher state contribution by the perturbative part 
above some threshold energy. For three-point Green’s functions, this assumption is more 
ambiguous than that for the two-point case because there are two energy variables € and č. 
It was argued by Blok and Shifman in [80], that the perturbative and the hadronic spectral 
densities cannot be locally dual to each other. The necessary way to restore duality is to 
integrate the spectral densities over the “off-diagonal” variable €_ = (€ — &')/2, keeping the 
“diagonal” variable €, = (€ + &)/2 fixed. It is w.r.t. €, that the quark-hadron duality is 
assumed for the integrated spectral densities. With this procedure the sum rule for the IW 
function yields [I6], 


= 4 1 a (qq)? "à (Q4GG), o 2w +1 
9f2pe- NT — _*_ i d & /T — B (w+1) s T? 10 
f Ee merh & p e qo cpm UE: (orig (10) 


where $ is the threshold energy for the sum rule. It should be noted that the Borel 
parameter has been chosen such that €(w = 1) = 1. The w dependence in the gluon 
condensate term and the coefficient of the exponential in the quark condensate term are 
different from Ref. [I5], but are consistent with Refs. [I7] [20]. From Eq. (10) one obtains 


the sum rule for the slope parameter of the IW function, i.e. p? = -gg =t It is given by 
w 
r 3 fË esar moa -™ | (0,GG) 
2 2 —A/T _ 5,—ec/T 0 E 2 
ape e desee e erac S. — m 


The baryonic decay constant f can be obtained by the sum rule Eq. (5). Note that we 
have not included the perturbative O(a,) corrections in the sum rule which is expected to 
be largely cancelled in the ratio of Eqs. (I0) and (5), or Eqs. (11) and (). This is what 
happened in the case of heavy meson form factors [8]. 

For the numerical analysis of the sum rule (II) we use the two-point sum rule () to 
eliminate the explicit dependence of Eq. (I) on f and A. This procedure reduces the 
uncertainties in the calculation. For the condensate contributions we take the standard 


values 
(gq) = —(0.24 + 0.01)? GeV? , (aa) = (0.012 + 0.004) GeV* (12) 


and mg = 0.8 GeV?. The threshold energies are taken to be equal, i.e. s = § = se. Imposing 
the usual criterion on the ratio of contributions of the higher-order power corrections 
and that of the continuum and using the central values of the condensates given in (12) 


there is an acceptable window of stability in the range 7' — 0.4 — 0.7 GeV in which the 
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Figure 1: The slope parameter of the IW function p? versus the Borel variable T. The 
curves refer to choices of the threshold parameters: se = 2.1, se = 1.9, se = 1.7, from top 
to bottom. 


calculated results do not change appreciably if the threshold parameter s, lies in the range 
1.7 < s, « 2.1 GeV. In Fig. 1, the sum rule for the slope parameter of the IW function 
p? is plotted as a function of the Borel parameter T for various choices of the continuum 
threshold in the range 1.7 < se < 2.1 GeV. One can see that the variation is very moderate 
for the Borel parameter in the range 0.4 « T' « 0.7 GeV. Our prediction for the slope 


parameter p? is given by 
p = 1.35 +0.12. (13) 


The errors reflect the uncertainty due to the sum rule window. The value of the slope 
parameter is in agreement with the recent experimental result [2] and the value obtained 
in a lattice determination [BI]. 

It is meaningful to ask why the present sum rule result is so different from those 
obtained previously. First, this maybe due to the systematic uncertainty of QCD sum 
rules. The systematic error resulting from the use of quark-hadron duality above s, is 
difficult to estimate. Conservatively speaking, there is a 1096 — 30% systematic error. 
Second, the previous linear fit to the IW function [16] [5] may be a rather poor fit in the 
w range: 1 — 1.43. In the following decay rate calculation, we assume an exponential form 


to parametrize the baryonic IW function, as the experiment did [2]: 


E(w) = &(1) exp[- e^ (v — 1)] . (14) 


Once we have computed the IW function we are now in the position to calculate the rate 


for the semileptonic A, — A, transitions. Neglecting the lepton mass the differential 


decay rate can be written as (see e.g. [4] [7]) 
dI u Gil Val mi r’ 

do /— 240? 
+o + 1) «Gt + (1 — 76: - w = D(62 GP) (15) 


vu — He —1) [sr +[(l+rFit w+ 1)(rFo t+ Fs)? 


where r = ma,/ma, and & = 1+ r? — 2rw. The form factors F; and G; can be expressed 
by a set of Isgur-Wise functions at each order in 1/mg in HQET. Taking into account the 
1/mg corrections, they are (see e.g. [A M) 


] A 
hate eee ee. Ee; E = E(w) 

-w-—l 25 À 
G =[1+(c+e)A—] lw), Fe=-Gs=-— Se), (10) 


where £g = 1/(2mg). Notice that the subleading Isgur-Wise function associated with the 
insertion of the Agcp/m, kinetic operator of the HQET Lagrangian has been neglected 
since it is negligibly small [I6] [I8]. In the numerical calculation we take the heavy quark 
masses to be my, = 4.8 GeV, m, = 1.4 GeV and A = 0.79 GeV Hj. With the masses of 
A, and A, given by the Particle Data Group (PDG) [82] the upper limit of w is wmar = 
(1 +1r?)/2r = 1.433. The decay rate can then be calculated to be 


T(A — Adr, = 242x10 !]|V4|? GeV . (17) 


The contribution of the 1/mg corrections amount to about 10%. 

Let us compare our theoretical prediction to experimental data. Recently, the DEL- 
PHI Collaboration measured the A, — Ar; branching ratio [2] as B(A, > A.€i%) = 
(5.0*li(stat)*t$(sys))06. On the other hand, the branching ratio given by the PDG is 
B(M, > ALl v; + anything) = (9.2 +2.1)% [82]. An error weighted average value was given 
by Albertus et al. [33]: 


The total A, decay width is determined by the inverse lifetime tx," where we take Ta, = 
1.229 + 0.80 ps [82]. The value of the CKM matrix element |V.4| can then be extracted to 
be 


[Va] = 0.041 + 0.004 , (19) 


where the error is due to experimental uncertainties. The value of [V.4| is in excellent 
agreement with the recent experimental determination by the DELPHI Collaboration 
[V4] = 0.0414 + 0.0012(stat) + 0.0021(syst) + 0.0018(theory), obtained from semileptonic 
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B — D*tr, decays [26]. It is also in good agreement with the value of |V.4| obtained from 
a nonrelativistic quark model calculation of the A, — A. transition [33]. 

In conclusion, we have presented a HQET sum rule analysis for the slope parameter p? 
of the baryonic IW function. We obtained p? = 1.35-£0.12 which is in good agreement with 
the recent DELPHI measurement and the results of a lattice calculation. When combined 
with the error weighted average value for the branching ratio of Ay — Ar; quoted in 
Eq.(I8), our integrated decay width including 1/mg corrections leads to a value for |V.| 
in excellent agreement with a recent determination of |Veæ| from B — D* decays. One 
should bear in mind that 1/ me corrections can decrease the decay rate by a few percent 
as estimated in [22]. They should be studied within HQET sum rules in future works. 
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